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Vortices in two-dimensional superconductors with broken time-reversal and spin-rotation symme- 
try can bind states at zero excitation energy. These socalled Majorana bound states transform a 
thermal insulator into a thermal metal and may be used to encode topologically protected qubits. 
We identify an alternative mechanism for the formation of Majorana bound states, akin to the way 
in which Shockley states are formed on metal surfaces: An electrostatic line defect can have a pair 
of Majorana bound states at the end points. The Shockley mechanism explains the appearance of a 
thermal metal in vortex-free lattice models of chiral p-wave superconductors and (unlike the vortex 
mechanism) is also operative in the topologically trivial phase. 

PACS numbers: 73.20.At, 73.20.Hb, 74.20.-z, 74.25.fc 



Two-dimensional superconductors with spin-polarized- 
triplet, p-wave pairing symmetry have the unusual prop- 
erty that vortices in the order parameter can bind a non- 
degenerate state with zero excitation energy [lHj]. Such 
a midgap state is called a Majorana bound state, because 
the corresponding quasiparticle excitation is a Majorana 
fermion — equal to its own antiparticle. A pair of spa- 
tially separated Majorana bound states encodes a qubit, 
in a way which is protected from any local source of deco- 
herence Since such a qubit might form the building 
block of a topological quantum computer there is 
an intensive search [7H12|| for two-dimensional supercon- 
ductors with the required combination of broken time- 
reversal and spin-rotation symmetries (symmetry class 

^ [H). 

The generic Bogoliubov-De Gennes Hamiltonian H 
of a chiral p-wave superconductor is only constrained 
by particle-hole symmetry, a x H*a x = —H. At low 
excitation energies E (to second order in momentum 
p = —ihd/dr) it has the form 

H = A(p x a x +p y a y ) + (U(r) +p 2 /2m)a z , (1) 

for a uniform (vortex- free) pair potential A. The elec- 
trostatic potential U (measured relative to the Fermi en- 
ergy) opens up a band gap in the excitation spectrum. 
At U — the superconductor has a topological phase 
transition (known as the thermal quantum Hall effect) 
between two localized phases, one with and one without 
chiral edge states [l3 - ll7| . 

Our key observation is that the Hamiltonian ((T|) on a 
lattice has Majorana bound states at the two end points 
of a linear electrostatic defect (consisting of a perturba- 
tion of U on a string of lattice sites). The mechanism 
for the production of these bound states goes back to 
Shockley [llj]: The band gap closes and then reopens 
upon formation of the defect, and as it reopens a pair 
of states splits off from the band edges to form localized 
states at the end points of the defect (see Fig. [1]). Such 
Shockley states appear in systems as varied as metals 
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FIG. 1: Emergence of a pair of zero-energy MS states as the 
defect potential Uo + SU is made more and more negative, at 
fixed positive background potential Uo = 0.3. (All energies 
are in units of 7 = HA/ a.) The energy levels are the eigen- 
values of the Hamiltonian JT]) on a square lattice (dimension 
100 a x 100 a, j3 = ti 2 /2ma 2 — 0.4 7, periodic boundary condi- 
tions). The line defect has length 50 a. The dense spectrum 
at top and bottom consists of bulk states. 



and narrow-band semiconductors [19( , carbon nanotubes 
[2fj| . and photonic crystals [2l|. In these systems they 
arc unprotected and can be pushed out of the band gap 
by local perturbations. In a superconductor, in contrast, 
particle-hole symmetry requires the spectrum to be ±E 
symmetric, so an isolated bound state is constrained to 
lie at E = and cannot be removed by any local pertur- 
bation. 

We propose the name Majorana-Shockley (MS) bound 
state for this special type of topologically protected 
Shockley states. Similar states have been studied in the 
context of lattice gauge theory by Creutz and Horvath 
[22l . |23| , for an altogether different purpose (as a way to 
restore chiral symmetry in the Wilson fermion model of 
QCD 0). 

Consider a square lattice (lattice constant a), at uni- 
form potential Uq. The Hamiltonian fl} on the lattice 
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FIG. 2: Main plot: Closing and reopening of the excitation 
gap at Uo = 0.3, ft — 0.4 (in units of 7), for states with 
k x = (black solid curve) and k x = n/a (black dashed curve). 
The MS states exist for defect potentials in between two gap- 
closings, indicated as a function of Uo by the shaded regions 
in the inset. (The red solid and blue dashed curves show, 
respectively Uo + SUo and Uo + SU-^. The label T indicates 
the topologically trivial phase.) 




Uo + 5U 



FIG. 3: Closing and reopening of the excitation gap at 
Uo = —0.3, P — 0.4 (in units of 7), for states with k x = (red 
curves) and k x = n/a (black curves). The results were ob- 
tained from numerical calculations using a constant isotropic 
pair potential A (solid lines) as in Fig. [2] as well as a spatially 
dependent, anisotropic pair potential (A x (r), A y (r)) deter- 
mined self-consistently from the gap equation (dashed lines) 



has dispersion relation 

E 2 = [U + 2/3(2 - cos ak x - cos ak y )] 2 

+ 7 2 sin 2 ak x + 7 2 sin 2 ak y . (2) 

(We have defined the energy scales /3 = H 2 /2ma 2 , 7 = 
HA /a.) The spectrum becomes gapless for Uq = 0, —4/3, 
and —8/3, signaling a topological phase transition [25l ]. 
The number of edge states is zero for Uo > and Uo < 
—8/3, while it is unity otherwise (with a reversal of the 
direction of propagation at Uq = —4/3). The topologically 
nontrivial regime is therefore reached for Uq negative but 
larger than —8/3. 

We now introduce the electrostatic line defect by 
changing the potential to Uq + SU on the N lattice points 
at r = (na, 0), n = 1, 2, . . . N. In Figs. [1] and [2] we show 
the closing and reopening of the band gap as the defect 
is introduced, accompanied by the emergence of a pair 
of states at zero energy. The eigenstates for which the 
gap closes and reopens have wave vector k x parallel to 
the line defect equal to either or ±7r/a (in the limit 
N — > 00 when k x is a good quantum number). 

We have calculated that the gap closing at k x = 
happens at a critical potential 8U = 6Uq given by [26| 

f- v /t/o({/ + 4/?)+7 2 for Uo > 0, 
SUo = < y/U (U + 4/3) + 7 2 for U < -4/3, (3) 

I no finite value otherwise. 

The critical potential SU^ for closing of the gap at 
k x = ±ir/a is obtained from Eq. ([3]) by the replacement 
of Uo with Uo + 4/3. The MS states appear for defect po- 
tentials Uo + SU in between two subsequent gap closings, 
as indicated in the inset of Fig. [5J 



We conclude that MS states exist for any value of Uo- 
In contrast, Majorana bound states in vortices exist only 
in the topologically nontrivial regime 0, H3] ■ The index 
theorem [28] for the production of zero-energy modes by 
the vortex mechanism, which requires the topologically 
nontrivial phase, is therefore not applicable to the Shock- 
ley mechanism. 

Our reasoning so far has relied on the assumption of 
a constant pair potential A, unperturbed by the de- 
fect. In order to demonstrate the robustness of the 
Majorana-Shockley mechanism, we have performed nu- 
merical calculations that determine the pair potential 
self-consistently by means of the gap equation [26|, l29| . 
In Fig. [3] we show a comparison of the closing and re- 
opening of the band gap as obtained from calculations 
with and without self-consistency, in the relevant weak 
pairing regime (Uo < 0). The self-consistency does not 
change the qualitative behavior. In particular, the gap 
only closes at k x = n/a for the parameters chosen (c.f. in- 
set in Fig. [2]) and the self-consistent determination of A 
only shifts the critical potential SU slightly. 

In Fig. H] we demonstrate that the MS states are lo- 
calized at the end points of the line defect. The expo- 
nentially small, but nonzero overlap of the pair of states 
displaces their energy from to ±E (with corresponding 
eigenstates ip— = a x ip* + related by particle-hole symme- 
try). The unpaired Majorana bound states tpi and ip2 
are given by the linear combinations 

^i = |(l -i)i) + + \(l + i)ip-, (4a) 
4>2 = \(l + «)V>+ + |(1 (4b) 

shown also in Fig. [H These states are particle-hole sym- 
metric, ip\ t 2 = & x ipl 2 ) so the quasiparticle in such a state 
is indeed equal to its own antiparticle (hence, it is a Ma- 



3 



2|^+| 2 |^i 12 




FIG. 4: Probability density of the paired (tp+) and unpaired 
{ipi,tp2) Majorana bound states at the end points of a line 
defect of length 50 a, calculated for Uo ~ 0.1 7, Uo + SU = 
-I.37, P = 0.4 7 . 
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FIG. 5: Average density of states for a potential that fluc- 
tuates randomly from site to site (U = 0.01 7, AU — 27, 
/3 = O.27). The lattice has size 400 a x 400 a. The right inset 
shows the same data as in the main plot, over a larger energy 
range. The left inset has a logarithmic energy scale, to show 
the dependence p oc In \E\ expected for a thermal metal (red 
dashed line). 



jorana fcrmion). 

If the line defect has a width W which extends over 
several lattice sites, multiple gap closings and reopenings 
upon increasing the defect potential 
Uo + 5U = —(hkp) 2 /2m to more and more negative val- 
ues at fixed positive background potential Uq. In the 
continuum limit W/a — > 00, the gap closes when [261 ] 
qW = mr + v, n = 0, 1, 2, . . ., with q = [k 2 F - (mA) 2 ] 1 / 2 
the real part of the transverse wave vector and v e (0, 7r) 
a phase shift that depends weakly on the potential. (Sim- 
ilar oscillatory coupling energies of zero-modes have been 
found in Refs. [3(J[3l|.) The MS states at the two ends 
of the line defect alternatingly appear and disappear at 
each subsequent gap closing. 

So far we constructed MS states for a linear electro- 
static defect. More generally, we expect a randomly vary- 



ing electrostatic potential to create a random arrange- 
ment of MS states. To test this, we pick U(r) at each lat- 
tice point uniformly from the interval (U — AC/, U + AU) 
and calculate the average density of states p(E). The 
result in Fig. [5] shows the expected peak at E = 0. This 
peak is characteristic of a thermal metal, studied pre- 
viously in models where the Majorana bound states are 
due to vortices [32h34I | . The theory of a thermal metal 
[IH predicts a logarithmic profile, p(E) oc In \E\, for the 
peak in the density of states, which is consistent with our 
data. 

Without Majorana bound states, the chiral p-wave su- 
perconductor would be in the thermal insulator phase, 
with an exponentially small thermal conductivity at any 
nonzero U [3, 32, 35, 36]. Our findings imply that electro- 
static disorder can convert the thermal insulator into a 
thermal metal, thereby destroying the thermal quantum 
Hall effect. Numerical results for this insulator-metal 
transition will be reported elsewhere (37j . 

These results are all for a specific model of a chiral p- 
wave superconductor. We will now argue that our find- 
ings are generic for symmetry class D (along the lines 
of a similar analysis of solitons in a polymer chain JUj ) . 
Let p be the momentum along the line defect and a a pa- 
rameter that controls the strength of the defect. Assume 
that the gap closes at a — c*o and at p — 0. (Because of 
particle-hole symmetry the gap can only close at p = 
or p = zszhn/a and these two cases are equivalent.) For 
a near ao and p near the Hamiltonian in the basis of 
left-movers and right-movers has the generic form 



H{a) = 



(vo - 



vi)p 



-i(a 
-(vo - 



- «o) 
vt)p 



(5) 



with velocities < v\ < vq- No other terms to first order 
in p — —ihd/dx and a — ao are allowed by particle-hole 
symmetry, H(a) = —H*(a). 

The line defect is initially formed by letting a depend 
on 1 on a scale much larger than the lattice constant. 
We set one end of the defect at x = and increase 
a from a(—oo) < ao to a(+oo) > ao. Integration of 
H[a(x)]ip(x) — then gives the wave function of a zero- 
energy state bound to this end point, 



ip(x) 



y /vo/vi - 1 
\Zvo/vx + 1 



exp 



a{x') — ao 
\/vl - v\ 



dx' 



(6) 

This is one of the two MS states, the second being at 
the other end of the line defect. We may now relax the 
assumption of a slowly varying a(x), since a pair of un- 
coupled zero-energy states cannot disappear without vi- 
olating particle-hole symmetry. 

We conclude with an outlook. We have identified a 
purely electrostatic mechanism for the creation of Majo- 
rana bound states in chiral p-wave superconductors. The 
zero-energy (mid-gap) states appear in much the same 
way as Shockley states in non-superconducting materi- 
als, but now protected from any local perturbation by 
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particle-hole symmetry. An experimentally relevant con- 
sequence of our findings is that the thermal quantum Hall 
effect is destroyed by electrostatic disorder (in marked 
contrast to the electrical quantum Hall effect). A recent 
proposal to realize Wilson fermions in optical lattices [39| 
also opens the possibility to observe Majorana-Shockley 
states using cold atoms. 

Our analysis is based on a generic model of a 
two-dimensional class-D superconductor (broken time- 
reversal and spin- rotation symmetry). An interest- 
ing direction for future research is to explore whether 
Majorana-Shockley bound states exist as well in the other 
symmetry classes [13|. Since an electrostatic defect pre- 
serves time-reversal symmetry, we expect the Majorana- 
Shockley mechanism to be effective also in class Dili 



(when only spin-rotation symmetry is broken). That 
class includes proximity-induced s-wave superconductiv- 
ity at the surface of a topological insulator [4Cj| and other 
exp erimentally relevant topological superconductors (4ll — 

sl- 
it would also be interesting to investigate the braiding 
of two electrostatic defect lines, in order to see whether 
one obtains the same non-Abelian statistics as for the 
braiding of vortices Lll . 
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Deutscher Akademischer Austausch Dienst DAAD, by 
the Dutch Science Foundation NWO/FOM, by an ERC 
Advanced Investigator Grant, and by the EU network 
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Appendix A: Line defect in lattice fermion models 

We calculate the closing and reopening of the excita- 
tion gap upon introduction of a line defect in a lattice 
fermion model with particle-hole symmetry. First we 
treat the Wilson fermion model [24J considered in the 
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main text, and introduced in the context of topological 
insulators in Refs. OSEj]. Then, in order to demonstrate 
the generic nature of the results, we consider an alter- 
native lattice model, the staggered fermion (or Kogut- 
Susskind) model Ew4al . introduced in the context of 
graphene in Refs. [37L |49|. 



Wilson fermions 



The Wilson fermion model has Hamiltonian 



H 



Tl n,m (nearest ncighb.) 



n, >n 



(Al) 



Each site n on a two-dimensional square lattice (lat- 
tice constant a) has electron and hole states |e) and \h). 
Fermion annihilation operators for these two states are 
collected in a vector c„ = (c„, e , c nt h)- States on the same 
site are coupled by the 2x2 potential matrix £ n and 
states on adjacent sites by the 2x2 hopping matrix T nm , 
defined by 0, SI 



c _ {U n \ _ / P 7 e 

Cn ~ 1 -U n ) ' lnm - U e - W — -p 



(A2) 



Here U n is the electrostatic potential on site n and 
Onm G [0, tt] is the angle between the vector 
and the positive y-axis (so 6 mn = tt — 9 nm ). In the con- 
tinuum limit a — > 0, the tight-binding Hamiltonian (|A1|) 
is equivalent to the chiral p-wave Hamiltonian (TTJ, with 
P = h 2 /2ma 2 and 7 = HA/ a. 

It is convenient to transform from position to momen- 
tum representation. For that purpose we take periodic 
boundary conditions in the y-direction, so that the trans- 
verse wavevector (in units of 1 /a) has the discrete values 
h = 2ttI/N, l = -(N- l)/2, . . . , — 1, 0,1, .. . ,(N — l)/2 
(for an odd number N of sites in the y-direction) . The 
Fourier transformation from position to momentum rep- 
resentation is carried out by the unitary matrix with el- 
ements [T] n i = N~ 1 / 2 e >nkl . We take an infinitely long 
system in the x-direction, so the longitudinal wavevector 
k varies continuously in the interval (— tt, tt]. 

For a uniform potential, U n = Uq for all n, the Fourier 
transformed Hamiltonian Ho(k) has matrix elements 

[H (k)] w - 6 w £i(k), (A3) 
£i{k) = Uo&z + 2P<t z (2 — cosfc — cos A;;) 

+ 7(c x sin k + a y sin ki). (A4) 

The corresponding dispersion relation is 

E(k, hf = [Uq + 2/3(2 - cos k - cos h)] 2 

+ 7 2 (sin 2 fc + sin 2 fc/), (A5) 

cf. Eq. ©. 

A line defect at row no (parallel to the a;-axis) adds to 
Hq the perturbation 



[8H]w = N^e^^'-^SUa, 



The determinantal equation Det (Ho + 5H — E) = for 
eigenenergy E reads 

Det (1 + Fl5Ua z Fo{Ho - E)- 1 ) = 0, (A7) 

in terms of an 1 x N matrix J-q with elements [.Fo]i/ — 
2V-VV n ° fe '. Sylvester's theorem, Det(l+AB) = Det(l+ 
BA), allows us to rewrite the determinant in the form 

Vet{l + 8Ua z To{Ho- E)- 1 ^) = 0, (A8) 



which reduces to 



= Dot 1 + SUa, 



N 4^ 



Si(k)-E 



( 1 £i(k)+E \ 

Det (i + SU* z -Y: e ^ 2 _ e2 ) 



(A9) 



A zero-mode is a pair of states (one left-mover and 
one right-mover) at energy E = 0. This can only occur 
at k — or k — tt (because for any eigenenergy E at 
k there must also be an eigenenergy — E at —k). From 
Eqs. (|A4[) and (|A9[) we obtain the condition for such a 
zero-mode, 



1 



£/ + 2/3(1 + <5- cos fc;) 



N^f[U + 2/3(1 



1 

«7' 



- S — cos ki)] 2 + 7 2 sin 2 ki 

(A10) 

where 6 = if k = and 6 — 2 if k — tt. In the 
limit N — > 00 we may replace the sum by an integral, 
N- x J2i -> (^r 1 S-ndh, which can be evaluated by 
contour integration. The resulting critical value of SU is 
given in the main text [Eq. ([3]) and following] . 



2. Staggered fermions 

The staggered fermion model is a discretization of the 
Hamiltonian (TT|) without the p 2 term. It is formulated in 
Refs. S^-S^l in terms of the transfer matrix Ai m , which 
relates the transverse wave functions ^ m +i = ■Mm^m 
at columns m and m + 1 (parallel to the y-axis). For 
a line defect along the ir-axis, the transfer matrix is m- 
independent, so we can omit the column number m. 

The transfer matrix (at energy E) has the form 



M = 



1 -iX 



(All) 



1 + iX' 

X = {ijy 1 ^^ + \Eg x J - \iGyU). (A12) 

In reference to Eq. ([1]), the parameter 7 = hA/a for 
lattice constant a. The N x N matrices J and JC have 
nonzero elements 



(A6) 



*Jn,n 1) 



r — _ 1 



(A13) 
(A14) 
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FIG. 6: Main plot: Closing and reopening of the excitation 
gap in the staggered fermion model. The MS states exist 
for defect potentials in the shaded regions in the inset. (All 
energies are in units of 7.) 



while the potential matrix IA (for a line defect at row no) 
is given by 

)• (A15) 

In momentum representation, the matrix X has ele- 
ments 



X lv =A l S w -i(SU/2 1 )a y - n: ^- 



(A16) 



'4cos 2 (jfcj/2)' 

where we have defined 

Ai = ia z tan(fc;/2) + {E/2 1 )a x - %(JJ Q /2rf)a v , (A17) 

vi = N- 1/2 e inokl (l + e ikl ). (A18) 

The dispersion relation of the staggered fermions is 
tan 2 (fc/2) = A(k, fc ; ) 2 , with 

A(k, hf = (£/2 7 ) 2 - tan 2 (fc ; /2) - ((7 /2 7 ) 2 . (A19) 

An eigenstate at energy E and longitudinal wavevec- 
tor k is an eigenstate of X with eigenvalue — tan(/c/2). 
The determinantal equation Det[A" + tan(fc/2)] = can 
again be simplified using Sylvester's theorem. The result, 
analogous to Eq. (|A9|) . is 



SU . 


1 


E 

1 


1 


27 v 


N 


Ai + tan(fc/2) J 


su . 


1 


E 


Ai - tan(fc/2) 


27 " 


N 


A(k,h) 2 -tan 2 (fc/2) 



(A20) 

Because of the pole in the dispersion relation at k = 7r, 
the zero-mode now exists only at k = 0. The condition 
for this zero- mode, analogous to Eq. (|A10|1 . is 



N ^ 



U /2j 



(C/ /2 7 ) 2 +tan 2 (fc J /2) 



27 
5U' 



(A21) 



For — > 00 we may again transform the sum into an 
integral, and thus obtain the critical potential 



SU = 




2 7 if Uq > 0, 
2 7 if U < 0. 



(A22) 



Upon varying the potential Uo+5U of the line defect, at 
fixed bulk potential Uo, the closing and reopening of the 
gap thus happens at Uo + SU = — 2 7 sign (Uq) (see Fig.[S]). 
The inset shows the region in parameter space where the 
Majorana-Shocklcy states exist in the staggered fermion 
model. This phase diagram is much simpler than the 
corresponding phase diagram for Wilson fermions (Fig. 
[51 inset), because of the absence of the extra parameter 
P (which quantifies the strength of the p 2 term in the 
Wilson fermion model). 



Appendix B: Self-consistent determination of the 
pair potential 

In order to determine the pair potential self- 
consistcntly in a spatially non-homogeneous situation, it 
is necessary to allow for a position-dependent, anisotropic 
pair potential A(r) = (A x (r), A y (r)). The Hamiltonian 
then reads [29| 

H =| {A x (r),p x } a x + \ {A y (r),p y } a y 

+ (U(r)+p 2 /2m)a z , (Bl) 

where {•,•} denotes the anticommutator. In the dis- 
cretization of this Hamiltonian on a square lattice, the 
spatial dependence of A(r) is taken into account in the 
hopping between neighbors as an average value of A(r) 
on the two lattice points. 

When the pair potential is homogeneous, the lattice 
Hamiltonian has the spectrum 

E 2 = [U + 2/3(2 - cosak x - cosak y )} 2 

+ 7 2 sin 2 ak x + j 2 sin 2 ak y (B2) 

with j x = HA x /a, jy — hA y /a and j3 = h 2 /2ma 2 . 

The Hamiltonian must be solved self-consistently to- 
gether with the equation for the pair potential. These 
read [29| (with derivatives discretized on the lattice) 

Ixir) = -ig ^2 ( u n(x + a,y) - u n (x - a,y))v*(x,y) 

E n >0 

- Un(x, y) (v*(x + a,y) - v*(x - a, y)), 
7y0) =9 E ( u n(x,y + a) -u n {x,y + a))v* n {x,y) 



E n >0 

- u n (x,y) (v*(x,y + a) - v* (x, y + a)). 



(B3) 



Here u n and v n are the electron and hole component 
of the wave function, respectively, and assumed to be 
from the tight-binding model, i.e. they are dimensionless 
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and represent the probability amplitude per lattice point 

The coupling constant g must be chosen such that it 
gives the correct pair potential 7 in the bulk. It can be 
calculated as 



A zero-mode ip = (it, v) is a (doubly degenerate) eigen- 



7 1 r 

- =-2 / d(ak x ) 

9 T 



d(ak y ) sin(afc x ) u(k)v*(k) 



d(ak x ) j d(aky) sin(aky) u(k)v* (k) , 

(B4) 



where it(fe) and v(k) are the electron and hole coefficients 
of the plane wave solutions of the bulk lattice Hamilto- 
nian with E > 0. 

In the particular case of a system that is translationally 
invariant in x-direction, as is the case for an infinitely 
extended line defect, the gap equations can be written 
as: 

lx(r) = jjr ^ u n {k x ,y)v* n {k x ,y) su\{ak x ) 



E n >0,k x 



_9_ 

N x 



£ 

E n >0,k x 

(u n (k x , y + a)- U n (k x ,y + a)) v* n (k Xl y) 

- u n (k x , y) {v*(k x , y + a) - v* n {k x ,y + a)) 

(B5) 

summing over N x longitudinal momenta k x , and solving 
the tight-binding problem for each k x individually. 

The self-consistent solution of the tight-binding Hamil- 
tonian and the gap equation (|B5I) is obtained in an itera- 
tive procedure. In the iteration, we neglect the influence 
of the vector potential arising from local currents [29( as 
those effects are expected to be minor for the examples 
considered in this work. Furthermore, we also avoid ad- 
justing the chemical potential Uq to obtain a fixed num- 
ber of electrons in the system and instead use a large unit 
cell so that the bulk value of A is recovered away from 
the defect. 



Appendix C: Line defect in the continuum limit 

We calculate the closing and reopening of the exci- 
tation gap upon introduction of a line defect in the 
Hamiltonian (p}, which is the continuum limit (a — » 0) 
of the Wilson fermion lattice model of App. IA II The 
mode matching calculation presented here is the one- 
dimensional version of the two-dimensional calculation 

in Refs. turn EH- 

The line defect, of width W, is formed by the electro- 
static potential profile 



U(r) 




if \y\ > W/2, 
if \y\ < W/2. 



(CI) 



state of the Hamiltonian ([T]) at E = 0, p x 
zero-mode should thus satisfy 

(U +p y /2m)u — iAp y v, 
(U +Py/2m)v = iApyU. 

For uniform U the solution is a plane wave, 



tfW = e tk °°' v 



s, s 



±1, 



with transverse wave vector 



= (m/h)(isA + s'^-A 2 - 2U/m). 



0. The 

(C2a) 
(C2b) 

(C3) 

(C4) 



In the region \y\ < W/2 the zero-mode ^ is a superposi- 
tion of the four states , ip^ , rp 1-, ip For y > W/2 

two decaying states with Im k ss i > appear in the super- 
position, while for y < ~W/2 the other two states with 
Im k SS ' < appear. In total ip has eight unknown coef- 
ficients, which we determine by demanding continuity of 
ij) and dip/dy at y = W/2 and y = ~W/2. The deter- 
minant of this set of equations should vanish, in order to 
have a nontrivial solution. There is only a zero-mode for 
Uq > 0, U Q + 8U < -mA 2 /2, determined by 



tan qW 



2qqg 
I 2 ~ Qo ' 



(C5) 



We have defined 



q = (m/h) v/-A 2 -(2/m)([/ + «/), (C6) 
qo = {m/h)^A 2 + 2U /m. (C7) 

The MS states exist in between subsequent gap closings, 
as indicated in Fig. [7] (shaded regions). 



U [)«A 2 /2] 



2 4 6 8 10 12 14 




FIG. 7: The red solid curves are the solution of Eq. (|C5|) for 
W = 4S/mA. The MS states exist in the shaded regions. 
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